"Breathing" rogue wave observed in numerical experiment 
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Numerical simulations of the recently derived fully nonlinear equations of motion for weakly three- 
dimensional water waves [V.P. Ruban, Phys. Rev. E 71, 055303(R) (2005)] with quasi-random 
initial conditions are reported, which show the spontaneous formation of a single extreme wave on 
the deep water. This rogue wave behaves in an oscillating manner and exists for a relatively long 
time (many wave periods) without significant change of its maximal amplitude. 

PACS numbers: 47.15.km, 47.35.Bb, 47.ll.-j 



I. INTRODUCTION 



The rogue waves (rare extreme events on a relatively 
calm sea surface, alternatively called freak, killer, or gi- 
ant waves), for a long time a part of marine folklore, since 
1970's have been methodically documented by oceanog- 
raphers (see review Q for examples and some relevant 
statistics). From the viewpoint of nonlinear science, a 
rogue wave is an extremely nonlinear object — typically, 
the amplitude of a freak wave in maximum is close to 
the amplitude of the corresponding limiting Stokes wave, 
that is h/X w 0.10 . . . 0.14, where h is the peak-to-trough 
height, and A is the length of the wave 0. Therefore, for 
adequate quantitative investigation, this phenomenon re- 
quires fully nonlinear equations and accurate numerical 
methods. For two-dimensional (2D) potential flows with 
a free boundary, a very efficient numerical scheme has 
been developed recently by Zakharov and co-workers . 
The scheme is based on exact (l-l-l)-dimensional equa- 
tions of motion written for the surface shape and for the 
boundary value of the velocity potential in terms of the 
so called conformal variables (the corresponding exact 
2D theory is described in Rcfs. 4, 5, 6., 1,1, ,9, 10] ). The 
method extensively uses algorithms of the discrete fast 
Fourier transform (FFT). With applying this method, 
impressive computations have been performed, where a 
numerical giant wave developed due to the Benjamin- 
Feir (modulational) instability [ill IT^ from a slightly 
perturbed Stokes wave. The spatial resolution in these 
numerical experiments was up to 2 • 10^ points As 
to three-dimensional (3D) flows, unfortunately, a similar 
exact and compact (2-}-l)-dimensional form of equations 
is absent. Therefore "exact" 3D simulations arc currently 
based on the rather expensive boundary element method 
(BEM) and its modifications (see El El E IlS: and 
references therein). Since the underlying algorithms of 
BEM are quite complicated, the best practically achieved 
spatial resolutions on the free surface for essentially 3D 
waves are typicall y fe w tens multiplied by few tens, as in 
the recent works |3 • Definitely, this is not suf- 
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ficient to simulate large wave systems with dozens and 
hundreds waves, as it is necessary for practical applica- 
tions. We exclude here the approximate equations de- 
scribing wave envelopes |l7l |18j 29_, 2Q\ , because they are 
not appropriate in real situations when many random 
waves with very different wave vectors and amplitudes 
are excited. Other approximate equations, for instance 
the weakly nonlinear Zakharov equations pH (23. . 
are only good if the wave steepness is small, but this 
is clearly not the case for the extreme waves. However, 
though rogue waves are strongly nonlinear, and the wave 
steepness cannot serve anymore as a small parameter of 
the theory, nevertheless another small parameter may ex- 
ist in the system. Namely, practically important is the 
situation when relatively long (much longer than a typ- 
ical wave length) wave crests are oriented along a defi- 
nite horizontal direction. For such weakly 3D fiows, the 
advantages of the conformal variables are almost recov- 
ered, as it has been explained in Refs. [2H2^. In the 
cited papers, the noncanonical Hamiltonian description 
in terms of the conformal variables has been generalized 
from 2D to 3D potential inviscid flows with a free sur- 
face, and the asymptotic expansion of the Hamiltonian 
functional on the small parameter e = [Ix/lq)^ has been 
suggested, where Ix is a typical wave length, and Iq is a 
large transversal scale along the wave crests. In partic- 
ular, the first-order 3D corrections have been calculated 
explicitly. What is important, all linear operators com- 
ing into the equations are diagonal in the Fourier repre- 
sentation. Therefore a relatively high spatial resolution 
(16384 X 256 in Ref.j2^) for the corresponding numeri- 
cal algorithm has been possible due to the large number 
of local operations that result from the Fourier diagonal- 
ization. In Ref. ,,26] some numerical results have been 
presented, for non-random initial conditions and typical 
dimensionless wave numbers about 20. 

In the present work another numerical experiment is 
reported, which is more close to reality. Main wave num- 
bers now are about 50, and the computations start with 
a quasi-random initial state (shown in Fig.^. Concern- 
ing efficiency of the numerical implementation, it should 
be noted that with the FFTW library '27], it takes less 
than 2 min to perform one step of the Rungc-Kutta-4 nu- 
merical integration on an Intel Pentium 4 CPU 3.60GHz 
with 2048M memory, for the maximal possible spatial 
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FIG. 1: (Color online). Map of the free surface at t = 0. 
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FIG. 2: Maximum wave height versus dimensionless time. 



resolution 16384 x 512. Here a giant wave formation has 
been observed as well, but contrary to the previous com- 
putations i2|| and l2a |. this freak wave is not breaking, 
but it exists for many wave periods without tendency 
towards increasing or decreasing its maximal amplitude 
(which in this case is distinctly less than the limiting 
Stokes wave amplitude, see Figs. I2l3f) . During the life 
time, the rogue wave behaves in an oscillating manner, 
with the highest crest being alternately ahead or behind 
of the deepest trough. Observation of such kind of be- 
havior is important for better understanding of the rogue 
wave phenomenon. 



II. EQUATIONS OF MOTION 

Here it is necessary to present the equations that were 
simulated. Their detailed derivation and discussion can 
be found in Refs. psi. l2^. We use Cartesian coordinates 
a;, q, y, with y axis up-directed. The symbol z denotes the 
complex combination: z = x + iy. For every value of q, 
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FIG. 3: Wave profiles at t = 50. 

at any time moment t, there exists an analytical function 
z(u + iv, q, t) which determines a conformal mapping of 
the lower half-plane of a complex variable w = u -\- iv 
on a region below the free surface. A shape of the free 
surface is given in a parametric form: 

Z = X + iY ^ z{u, q,t) = u+{i- H)Y{u, q, t). (1) 

The Hilbert operator H is diagonal in the Fourier repre- 
sentation: it multiplies the Fourier-harmonics 

Yk^{t) = jY{u,q,t)e~'''^-'""'dudq 

by [i sign k] , so that 

HY{u,q,t) = / [isign/c]rfe™(i)e*'="+™«dfcdm/(27r)2. 



(2) 

Thus, the first unknown function is Y{u, q, t). The second 
unknown function is the boundary value ip{u,q,t) of the 
velocity potential, 

il}{u,q,t)= I i/;fc„(t)e*'="+™«dfcdm/(27r)^ 



Correspondingly, we have two main equations of motion. 
They are written below in a Hamiltonian non-canonical 
form involving the variational derivatives {SIC/ Sip) and 
(SIC/SZ), where IC{ip, Z, Z} is the kinetic energy. The 
first equation is the so called kinematic condition on the 
free surface: 



Zt = iZuil + iH) 



{SK/Sj}) 



(3) 



The second equation is the dynamic condition (Bernoulli 
equation) : 



-gim Z — ipuH 



[SIC /Sip) 

\z„p 



Im 



((1 - iH) [2{5}C/5Z)Zu + {5}C/SiP)^u] 



-,(4) 
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FIG. 5: (Color online). Spectrum of the right-propagating 
waves ai t — 50. Here shown is logj^Q \Ykm + iPkm\, where 
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FIG. 6: Distribution of the surface elevation y{x,q) at t = 50 
(no averaging over the time is done). 



parameter e ~ [Ix/lqY' <§; 1), we have an approximate 
kinetic-energy functional in the form 



FIG. 4: (Color online). Top: map of the free surface at f = 50 
(7 min 30 sec). The rogue wave has coordinates a; « 3.7 km, 
q « 3.7 km. Middle: map at t = 60 (9 min 1 sec). The rogue 
wave is at a; « 4.3 km, q ~ 3.4 km. Bottom: map at f = 91 
(13 min 40 sec). The rogue wave is at s ~ 1.3 km, q ~ 2.8 
km, and a specific wave pattern behind of it is visible. 



where g is the gravitational acceleration. Equations Q 
and Q completely determine evolution of the system, 
provided the kinetic energy functional /C{'0, Z, Z} is ex- 
plicitly given. Unfortunately, in 3D there is no exact 
compact expression for IC{ip^ Z, Z}. However, for long- 
crested waves propagating mainly in the x direction (the 



K. K K. — J ipHipududq + (5) 

where the first term describes purely 2D flows, and weak 
3D corrections are introduced by the functional 

# = ^ y {Zu^q - Zq^u)G{Zu-^q- Zq^u)dudq 
+ ^ y" I [{Zu-^q - Zq^uf/Zu] E{Z~U) 

-{Z-u) E[{Z^^q - Zq^^Y/Z^] I du dciQ) 
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Here ^ = {1 + iH)ip, and the operators G and E are 
diagonal in the Fourier representation: 



G(fc,m) = 



E{k, to) 



'2i 



2\k\ 



\kl 



(7) 



(8) 



A difference between the above expression (0) and the 
unknown true water-wave fcinetic energy functional is 
of order e^, since G{k^O) = l/{ik) for positive fc, and 
E{k,0) = 1 (see Refs. ISl^^). Besides that, the linear 
dispersion relation resulting from /C is correct in the en- 
tire Fourier plane (it should be noted that in Ref. 
another approximate expression for JC was used, also re- 
sulting in the first-order accuracy on e and correct linear 
dispersion relation). Thus, we have (SlC/dijj) ~ (SJC/Sip) 
and {5K./5Z) « {6J- /6Z) in equations H3I4|) . with explicit 
expressions closing the system: 



2 Re 



(9) 



6 
4 

>• 

-2 
-4 
-6 



n— 
H 

q= 


14/16)*5 km 

13/16)*5 km 










q= 


11/16)*5 km 


























































J 



































3.1 3.15 3.2 



(a) 



3.25 3.3 3.35 

X [km] 



3.4 3.45 3.5 



6 
4 

>• 

-2 
-4 
-6 



q= 
q= 
^ q= 
q= 


13/16)*5km - 
12/16)*5km - 
11/16)*5 km 
10/16)*5km - 






























■/ // 














\ / 























(b) 



3.5 3.55 3.6 3.65 3.7 3.75 3.8 3.85 3.9 
X [km] 



Jz 



-(*, - Zq^u/Zu)E{Z ~ u) 



Zu dq G{Zu^q - Zq-^u) 



+ (*, - Zg'^jZu)E{Z - u) , (10) 



T77 = -7;'^,du 



(*, - Z,*JZ,0S(Z - u) 



-(*, - Z,^'^jZ,,)EiZ - u) 



16 



du[{^^-Zq^ulZuYE{z-u)] 



-E{^^-Z^^ulZufZu .(11) 



III. RESULTS OF THE NUIVIERICAL 
EXPERIIVIENT 
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Following the procedure described in Ref. |2y|, a nu- 
merical experiment has been performed, which is de- 
scribed below. A square 5x5 km in (u, (7)-plane with 
periodic boundary conditions was reduced to the stan- 
dard square 27r x 27r and discretized by x L points. 



FIG. 7: Rogue wave profiles at f = 45, t — 50, t = 55, and 
t = 60. 
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Thus, all the wave numbers k and m are integer. Di- 
mensionless time units imply g = 1. As an initial state, 
a superposition of quasi-randomly placed wave packets 
was taken, with 25 packets having wave vector (60,2), 
25 packets having wave vector (50,0), 16 packets with 
(40, -2), and 12 packets with (30, 1). Amphtudcs of the 
packets with fc = 50 were dominating. Thus, a typical 
wave length was 100 m, and a typical dimensionless wave 
period T = 2it/V50 ~ 1. The crest of the highest wave 
was initially less than 3 m above zero level. A map of the 
free surface at t = is shown in FigQ] It is clear from 
this figure that initially e ~ 0.01. 

The evolution of the system was computed with N — 
16384 and i = 256 to i = 40.0, until beginning of a 
rogue wave formation. After t = 40.0, the rogue wave 
was present in the system (see Fig. 0J), and during many 
wave periods its height in maximum was approximately 
7 m, as Fig. 13 shows. It resulted in widening of the wave 
spectrum (see Fig. |31 where e ~ m? /k"^ ~ 0.05), and 
therefore L = 512 was employed from t = 40.0 to t — 
60.0. Within this period, the total energy was decreased 
by 0.5% due to numerical errors. Finally, from t = 60.0 
and to the end of the experiment, L = 1024 was used 
to avoid progressive loss of the accuracy (the last stage 
has required computer with 3072M memory, and it took 
5 min per one step of integration). 

The presence of rogue wave strongly affects the prob- 
ability distribution function P{y) of the free surface el- 
evation y{x,q). Fig. shows that the distribution has 
a Gaussian core and "heavy" tails, which are not sym- 
metric - large positive y are more probable than large 
negative y. 

The most interesting observation of the present nu- 
merical experiment is that the freak wave can exist for a 
relatively long time without significant tendency towards 
breaking or disappearing. While "living" , the big wave 
does something similar to breathing, as shown in Fig. \7\ 
The rogue wave propagates along the free surface (with 
the typical group velocity, but there is also a displace- 
ment in g-direction) , and position of the highest crest is 



alternately ahead or behind of the deepest trough. Very 
roughly this behavior corresponds to a short wave enve- 
lope (with approximately one wave length inside) filled 
with a strongly nonlinear periodic Stokes-like wave. The 
time period of this "breathing" roughly equals to two 
wave periods, which property seems natural due to the 
fact that the group velocity of the gravitational waves 
is one half of the phase velocity. After 11 periods of 
"breathing" with the almost constant amplitude 7 m, the 
rogue wave gradually irradiates accumulated energy into 
a specific wave pattern visible in Fig. 21 at t = 91. This 
wave pattern nearly corresponds to the resonance condi- 
tion 

u;(ko) + VgrCko) ■ (k - ko) - Lo(k) = 0. 

where the wave vector kg — {ko^mo) characterizes the 
rogue wave, and w(k) = (gjkl)^/^ is the linear disper- 
sion relation. However, a more accurate explanation and 
an analytical study of the observed coherent nonlinear 
structure is a subject of future work. 



IV. SUMMARY 

Thus, the recently developed fully nonlinear theory for 
long-crested water waves, together with the correspond- 
ing FFT-based numerical method '26] are shown in this 
work to be an adequate tool for modeling rogue waves in 
close to real situations, that is with many random waves 
propagating mainly along a definite horizontal direction. 
Now it has been possible to deal with quite high spatial 
resolutions, since in the present algorithm all the non- 
local operations are reduced to the FFT computing, and 
the latter is really fast with modern numerical libraries. 
Different dynamical regimes of the rogue wave evolution 
can be investigated. In particular, the present article re- 
ports observation of a long-lived rogue wave. Such waves 
are definitely important from practical viewpoint. 
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